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ABSTRACT 

We investigate the sampling and dipole convergence properties of flux-limited 
samples of mock X-ray clusters in relation to their underlying "parent" cluster 
distribution. To this purpose, we resort to numerical simulations of the cluster 
distribution and extract samples resembling the main observational features of 
X-ray selected cluster samples. The flux-limited samples, being quite sparse, 
underestimate the amplitude of the "parent" cluster dipole by ~ 15 per cent 
on average for Local Group-like observers. However, the general shape of their 
dipole amplitude profiles are in relatively good agreement. We also calculate 
the expected contribution of clusters, selected according to the relevant crite- 
ria, to the soft (i.e. 0.1 — 2.4 keV) extragalactic X-ray background, using the 
ESO Key Project X-ray luminosity function, assuming a flat universe with 
vanishing cosmological constant. We obtain a value of about 10 per cent of 
the observed XRB flux. 

Keywords: X-ray clusters: clustering - large scale structure of Universe 
- dark matter - soft X-ray Background 



1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally- 
collapsed structures in the universe. This property 
makes these objects potentially useful tracers of the 
global matter distribution. Considerable observational 
effort has thus been spent in compiling cluster cata- 
logues as free as possible from errors and selection bi- 
ases. While it remains a very difficult task to construct 
bias-free optically-selected samples of rich clusters, se- 
lecting objects through their X-ray emission offers a 
cleaner possibility (e.g. Cen 1997). The powerful X-ray 
emission from the hot intra-cluster gas, peaked at the 
gravitational centre of the clusters potential wells, de- 
fines them much better than optical light and it is more 
directly related to their total mass, making such X-ray 
emitting objects prime candidates for tracers of the mat- 
ter distribution on very large scales (for recent reviews 
see Ebeling et al. 1995; Guzzo et al. 1995; Collins et al. 
1995; Gioia 1996). 

The analysis we shall undertake in this paper was 
initially inspired by the ROSAT all-sky survey (RASS) 
and the follow-up work forming the ESO Key Program 
on southern cluster redshifts (Guzzo et al. 1995; de 
Grandi 1996). A number of other X-ray cluster surveys 
based on ROSAT are also available; the South Galac- 



tic Pole project (SGP) by Romer et al. (1994) (a flux- 
limited sample comprising 128 X-ray clusters restricted 
to an area of 3100 deg 2 ), the extensive optical follow-up 
work aimed at selected clusters from RASS by Allen et 
al. (1992) and by Crawford et al. (1995) (that accom- 
panied and led to the X-ray compilation of the brand 
new Brightest Cluster Sample (BCS) by Ebeling et al. 
1997b) and last but not least at all, the recent X-ray 
Brightest Abell-type cluster sample (XBACs), which 
constitute the largest all-sky flux-limited cluster sam- 
ple to date (Ebeling et al. 1996b). Prior to the ROSAT 
mission many other X-ray projects were conducted with 
the aim to construct a large unbiased X-ray flux-limited 
cluster sample that could be used for detailed cosmolog- 
ical and statistical purposes (mapping the large scale 
structure of the universe, quantifying cluster evolution, 
investigating correlation function properties, establish- 
ing the X-ray luminosity function for clusters of galax- 
ies, to name but a few). 

The principal aim of this work is to investigate 
the disparities of dipole amplitude, structure and align- 
ment between the X-ray flux-limited samples and their 
"parent" cluster distribution. We shall investigate this 
using 'mock' cluster samples generated in numerical 
simulations via an optimized version of the Truncated 
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Zel'dovich approximation (TZA; Coles, Melott & Shan- 
darin 1993; Borgani, Coles & Moscardini 1994; Bor- 
gani et al. 1995). This study therefore complements 
and extends that of Kolokotronis et al. (1996) and of 
Tini Brunozzi et al. (1995) (TB95 hereafter) which con- 
cerned the relationship between the dipole of optically- 
identified clusters, that of galaxies, and that of the un- 
derlying matter distribution in simulations. The authors 
found that the high peaks of the density field (clusters) 
do trace the underlying mass dipole structure reason- 
ably well, with a roughly constant amplitude difference 
corresponding to a linear biasing term, though with a 
relatively large observer-to-observer scatter. The degree 
of cluster-mass dipole correlation depends on the den- 
sity threshold above which the clusters are defined and 
for the relatively low density threshold, which corre- 
sponds to a mean separation of ~ 30/i -1 Mpc (similar 
to that of the ESO KP X-ray clusters), they found a 
very good correlation. However, the above cluster-mass 
dipole correlation starts degrading with increasing in- 
tcrcluster separation, or density threshold (see section 
4.1 of Kolokotronis et al. 1996). 

Dipole analyses have been carried out by many 
authors in recent years; for a complete description of 
the literature, see the references in Kolokotronis et al. 
(1996) (hereafter K96). 

Using linear perturbation theory (Peebles 1980) one 
can relate the cluster dipole and the observer pecu- 
liar velocity to obtain an estimate of the /^-parameter 
(= n[)' 6 /6 c , where b c is the cluster bias parameter with 
respect to the mass). It is worth noting here specifically 
that, different classes of extragalactic objects (QSOs, 
AGN, galaxies, clusters of galaxies) have different rela- 
tionships to the underlying matter distribution and may 
therefore differ substantially in their respective dipole 
behavior. Although this complicates the use of dipole 
properties as estimators of the cosmological density pa- 
rameter, Qo, it does allow one in principle to study the 
relative bias displayed by such objects (K96; Plionis 
1995). Prominent among those analyses that have used 
X-ray selected objects as tracers of the density field are 
Miyaji & Boldt (1990) (X-ray AGN) and Lahav et al. 
(1989) (X-ray clusters). 

Our subsidiary intention, having generated the 
mock X-ray cluster samples complete down to a certain 
flux limit, is to estimate their contribution to the soft X— 
ray Background (XRB) which, ever since its discovery 
(Giacconi et al. 1962), has been one of the outstanding 
puzzles in this field. Such an application is particularly 
interesting because as far as we know we will be simulat- 
ing a large sample of X-ray clusters (for similar X-ray 
samples see also McKee et al. 1980; Piccinotti et al. 
1982; Maccacaro et al. 1984; Lahav et al. 1989; Gioia et 
al. 1990; Edge, Stewart, Fabian & Arnaud 1990; Edge 
& Stewart 1991; Henry et al. 1992; Guzzo et al. 1995; de 
Grandi 1996; Ebeling et al. 1995; Ebeling et al. 1997b 
and references therein). This enables us to place tight 
limits on the relation between X-ray emitting clusters 
and the XRB in the soft energy band (0.1-2.4 keV). We 
also stress the fact that, a substantial number of work is 
dedicated to measure the contribution of X-ray objects 
to the hard XRB (2-10 keV). Although the availability 



of data from the ROSAT satellite makes this task feasible 
for energies below 2 keV, the fact that both the spec- 
trum of the XRB and the contribution of the resolved 
discrete component are yet to be determined, forces us 
to make some relatively uncertain assumptions, as we 
point out later. 

The outline of the paper is as follows. An overview 
of the dipole analysis, together with techniques for an 
accurate determination of the inherent shot-noise errors, 
is presented in Section 2. We present a useful set of 
definitions, observational constraints and the method of 
reproducing reliable simulations of the relevant X-ray 
populations in section 3. We discuss the cluster dipole 
results in section 4, while an estimate of the contribution 
of X-ray clusters to the soft XRB in a critical density 
model is presented in section 5. Finally, we present our 
main conclusions in section 6. 



2 DIPOLE FORMALISM 

In this section we will present the method we use to 
calculate the peculiar gravitational acceleration induced 
by the distribution of matter as it is traced by clusters 
of galaxies. 



2.1 Theoretical Preamble 

As we discussed more extensively in K96, we need to as- 
sume the standard gravitational instability picture for 
all of this work. Within this picture we assume that lin- 
ear perturbation theory applies and, perhaps with less 
justification, that a linear biasing scheme relates fluctu- 
ations in cluster numbers to fluctuations in the overall 
density of material, i.e. 



(i) 



(Kaiser 1984; Coles 1993). This scheme motivates the 
introduction of the parameter /3, defined by 



/(n„ 



b ' 



(2) 



where / is the mass fluctuation growth rate in linear per- 
turbation theory (e.g. Peebles 1980). Using linear the- 
ory we can relate velocities and accelerations. A clas- 
sical case is to relate the only very accurately known 
observed peculiar velocity in cosmology, that of the Lo- 
cal group (LG) relative to the cosmic microwave back- 
ground (CMB), u(r), to the peculiar gravitational ac- 
celeration induced upon it by the surrounding density 
fluctuations, g(r), which then, provided that the two 
vectors are aligned, allows (3 to be estimated via: 



i [ 5 ^w dr > 



(3) 



where S(r) is the density fluctuations of the underlying 
mass distribution at position r. Knowledge of b allows 
an estimate of /3, if Q is given as an input parameter 
and vice versa. Although the assumption (1) is rela- 
tively well-motivated in a statistical sense, it may well 
turn out to be a poor approximation when applied on 
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a point-to-point basis, as required by this analysis. In 
fact, K96 found by comparing mass and cluster dipoles 
in simulations, in which fio was an input parameter, 
that they recovered the correct mean bias value but al- 
beit with a large scatter (see also section 4.4). 

Note that we will be using distance units in km s _1 
and thus the H a dependence drops. Wherever necessary, 
we will use H a = 100/ikms -1 Mpc -1 . 



2.2 Analysis Tools 

We employ the method of moments to quantify the dis- 
tribution of clusters around the observer. Having an 
available set of cluster positions, n, around an observer 
placed at the origin, we can calculate the dipole and 
monopolc moments as follows: 



N(R) 

m = wj 



and 

N(R) 

D — S Wifi 

i=l 



(4) 



(5) 



The sum is taken over all N(R) clusters within a radius 
R and Wj (oc rr 2 ) are the appropriate weights, which 
should also take into account any sample selection ef- 
fects. The monopole should increase linearly with dis- 
tance R while the dipole D keeps increasing until the 
most distant density inhomogeneity, that affects the dy- 
namics of the observer, is taken into account. The scale 
at which this happens, which we denote as the con- 
vergence depth (f?conv), should be within the effective 
depth of the catalogue if the estimated value of f3, via 
equation ^ is to be reliable. 

Since the samples we are interested in are flux- 
limited, the number of objects decreases with radial dis- 
tance due to the rapidly decreasing selection function 
4>(r), defined as the number density of objects that can 
be seen above the limiting luminosity at distance r: 



0(r) = 



$(L) dL, 



(6) 



where is the luminosity function of the objects 

under study and L m i n (r) = 47rr 2 5'ii m , with Sum being 
the flux limit of the sample under study. In such a case 
the simple gravitational weights should be replaced by 
wi oc r ; -2 , where the inverse selection function 

corrects (statistically) for the unseen portion of the lu- 
minosity function (for more details see K96 and refer- 
ences therein). 

A word of caution is due here. Since we sample clus- 
ters up to 0.16 in redshift, we should have taken into 
account possible evolutionary effects of the X-ray lumi- 
nosity function. Instead, we considered that the X-ray 
luminosity function (XLF hereafter) holds equally well 
for local and for distant objects which is clearly a sim- 
plification. Although this is still an open issue, there is 
recent evidence (Ebeling et al. 1996a; Nichol et al. 1997; 
Collins et al. 1997; Ebeling et al. 1997a; Romer et al. 
1997 and references therein) that the cluster XLF does 



not evolve at least up to z ~ 0.2-0.3. Therefore we feel 
that using an non-evolving XLF is a safe approximation 
at least for the distance range under study. Had we used 
an evolving XLF, this would obviously have had a direct 
impact on the accuracy of estimates of the contribution 
of X-ray clusters to the soft XRB (depending on the 
degree and the nature of the evolution) . 

The combination of the observed and predicted pe- 
culiar velocity yields an estimate of the /3-parameter via 
equation (^[) which is equivalent to: 



u(r) = —D 

47rn 



/3V C , (7) 

where D is estimated via equation (|^) and the sub- 
script c corresponds to cluster velocity. The importance 
of using an accurate estimator for the mean density of 
clusters is evident; there are alternative ways of con- 
structing such an estimator and we refer the reader to 
K96 for details (Section 2.1). We verified that using dif- 
ferent density estimators, the outcome is within 10 per 
cent of each other (see also Davis & Huchra 1982). In 
what follows we just estimate densities using equation 
([]). Note that this method only generates a reliable esti- 
mate of 13 if (i) the two vectors u and D are well-aligned 
(say within 25°) and (ii) that the effective depth of the 
sample, -Rmax, exceeds R CO nv However, we should bear 
in mind that the very local contributions to the dynam- 
ics of the observer may not be represented in the cluster 
distribution, since these can be due to nearby galaxies, 
groups of galaxies or even voids in the local matter dis- 
tribution. Such local effects can also be viewed as the 
cause of the observer-to-observer scatter in the linear 
biasing relation between the cluster and mass distribu- 
tions (TB95; K96) 

2.3 Shot Noise Effects 

The sparseness with which the flux-limited cluster sam- 
ples trace their underlying parent population introduce 
shot-noise (discreteness effects) in their dipole estimates 
which increase as a function of distance. The method of 
angular position reshuffling, used in K96 to estimate a 
similar effect for the IRAS galaxies, is not suitable in 
this case. Had we used the positional reshuffling method, 
we would have been measuring the convolution of the 
shot-noise introduced from the dilute sampling of the 
parent cluster population by its flux-limited subsample 
and of a sort of intrinsic parent cluster population shot- 
noise dipole (under the false assumption that clusters 
are Poisson samplers of the underlying mass field). In 
reality, since we will be inter-comparing cluster dipoles 
we are interested only in the former type of shot-noise 
which we estimate according to the formalism developed 
in Strauss et al. (1992): 



IDI; 



1-1 



(8) 



where the sum is over all clusters in the flux-limited 
sample. Evidently, the shot-noise dipole increases with 
depth and thus the amplitude of the flux-limited cluster 
dipole will also have a similar behavior, even beyond the 
true underlying dipole convergence depth. Therefore, in 
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order to identify this convergence depth, a parameter 
which is very important for a reliable determination of 
(3 via equation ) , we must correct the raw flux-limited 
dipole for the effects of shot-noise. Although the ampli- 
tude of the shot-noise dipole is easy to estimate, such a 
correction is not straightforward, because the shot-noise 
dipole direction is random and could point anywhere. 

In the limit of uniform sampling we have from the 
central limit theorem that the three Cartesian com- 
ponents of the shot-noise dipole are equal and thus 
.D|sn — \D\sn/V3- Using a coordinate system in which 
the ID shot-noise dipole is parallel to the direction of the 
shot-noise free dipole we have that |£>| cor r ~ |-Dj ra w ± 
|-D|sn- We will correct for the effects of shot-noise using 
|-D| 0OII ~ | -D | raw — \D\ln to provide maximum discrete- 
ness correction and thus a rough lower limit to the es- 
timated cluster dipole while we will also investigate the 
case of no shot-noise correction which should provide 
an upper limit to the estimated dipole. The classical 
shot-noise correction, |X? C or| = -\/|-Draw| 2 — |-D?n| 2 al- 
ways produces results intermediate between the above 
limits and therefore we will not consider it further. 

Finally, an estimate of the misalignment induced by 
the shot-noise dipole on the true underline dipole direc- 

3 ]ETr 



tion is given by: 88 sn w arctan 



(see also 



Juskiewicz, Vittorio & Wyse 1990 and Lahav, Kaiser 
Hoffman 1990 for dipole misalignment correction). 



3 SAMPLES AND SIMULATIONS 

We now turn to the observational properties of the X- 
ray cluster samples available at the present time, and 
the methods we use to simulate them. For more details 
on the observational data and motivation for this work, 
see Guzzo (1995), Guzzo et al. (1995) and de Grandi 
(1996). 

3.1 Properties of the Observed X— ray Cluster 
Samples 

In order that our simulations be as closely related as 
possible to the available X-ray cluster data, we need 
to use an appropriate X-ray luminosity function and a 
realistic flux-limit. We use the Schechter-like form for 
the luminosity function as provided by de Grandi (1996) 
with parameters given in Table 1: 



$(L) 



(9) 



Here, L„ is the characteristic luminosity, A is the overall 
normalization of the number-density, and a is the usual 
power-law index. This luminosity function has been 
found to describe adequately a pilot sample comprising 
111 X-ray clusters having redshifts between 0.02 and 0.2 
above the limiting flux of 3 x 10 -12 ergs -1 cm -2 . The 
energy band for this application is 0.5-2 keV, instead of 
the broader 0.1-2.4 keV, which we will be using. How- 
ever, we can easily convert results (fluxes or luminosi- 
ties) from the hard band to the band merely by using 
the conversion factor of order of 1.615, as suggested by 
de Grandi (private communication; see also section 3 of 



Table 1. Parameters for the de Grandi (1996) luminosity 
function 

Parameters Values Units 

A 2.3153 x 10 -6 (/i -1 Mpc) -3 (10 44 /i -2 ergs -1 ) 0-1 

1.32" 1 



a 



,+0.21 

-0.23 

0.66l°; 2 ^ x 10 44 



de Grandi 1996 for details). Note that the above param- 
eters of the luminosity function were derived from the 
original parameters of de Grandi (1996), by using the 
parametrised form of the Hubble constant. We prefer 
this conversion and consequently report results, which 
are /i-dependent (see sections 4 and 5). 

It is necessary to implement a minimum distance 
criterion, because nearby clusters are too extended to 
be identified by the present analysis software (SASS) of 
the ROSAT survey. Furthermore, the SASS detection al- 
gorithm implementation leads to an underestimation of 
the fluxes for all kinds of X-ray sources. In fact, both the 
pilot sample of the ESO KP and the whole ESO KP sur- 
vey make also use of an alternative flux estimator (for 
both pointlike and extended X-ray emission), namely 
the steepness ratio technique (for a detailed delineation, 
see de Grandi et al. 1997 and references therein). Such 
a method greatly improves the flux incompleteness ini- 
tially produced by the SASS, especially when the signal- 
to-noise ratio for the RASS sources is small (i.e. when 
photon statistics is low). Nevertheless, we adopt a value 
of iCin ^ 45/i -1 Mpc (z min « 0.015 as in RASS) as 
the minimum fiducial distance at which the sample is 
complete. Although the maximum luminosity used for 
the pilot sample is 1.01 x 10 45 ergs -1 , the most lumi- 
nous X-ray cluster identified within RASS yields the 
enormous value of 6.2 x 10 45 ergs -1 (Schindler et al. 
1995). This object appears to be incompatible with the 
luminosity function (equations 6 and 11). In any case, 
this is a preliminary estimate of this function, so that 
our results which are based on it, should be taken as 
indicative rather than definitive. Since the integration 
of this luminosity function is highly insensitive to L max 
for values well above L t , we impose an upper limit of 
1 — 1.2 x 10 45 h~ 2 ergs -1 . The lower luminosity limit can 
be easily found in conjunction with J? m i n and the RASS 
flux limit. Using the ESO Key Project flux limit, i.e. 
Slim = 3 x 10 -12 ergs -1 cm -2 , produces samples of few 
clusters, spanning a large redshift range, which there- 
fore provide unreliable dipole estimates, principally be- 
cause of the very large shot-noise uncertainties. Since 
we are interested in producing theoretical predictions 
that can be compared with future, more sensitive sur- 
veys we use a deeper flux-limit, comparable to that of 
the RASS SGP project, set by Romer et al. (1994) and 
those by Collins et al. (1995) and Guzzo et al. (1995) 
i.e. Si im = 1.5 x 10 -12 
to L min k4x 10 41 h 
hope that a survey with this deeper flux limit would still 
be adequately described by the same luminosity func- 
tion. Since, the pilot catalogue as given above, reflects a 
quite uniform picture of the whole sampling space, it is 
our reasonable hope that such an alteration of the flux 



erg s 1 cm A , which corresponds 
2 erg s -1 . We do, however, have to 
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limit will not affect the parameters of the luminosity 
function (L. Guzzo, private communication). 

3.2 The simulations 

We now briefly describe the simulations we use to gener- 
ate mock samples of galaxy clusters. For a more detailed 
description, see Borgani et al. (1995). 

Our prescription is based on an optimized version 
of the Truncated Zel'dovich approximation (Zel'dovich 
1970; Shandarin & Zel'dovich 1989; Coles, Melott & 
Shandarin 1993). In this approximation, known as TZA, 
particles move along straight lines in response to the 
initial gravitational potential generated by a smoothed 
version of the initial density field. This describes the dy- 
namics quite accurately up to the shell-crossing point, 
which is suppressed by the initial smoothing. In the 
simulations we use here, density and velocity poten- 
tials are reconstructed on 256 3 grid points in a box of 
L = 960ft -1 Mpc aside. It has been shown that TZA 
accurately locates clusters when as < 1 (Borgani et al 
1995; Moscardini et al. 1996). Indeed, this is the case for 
the model we will be dealing with in this paper, a cold + 
hot dark matter model (CHDM hereafter) with f2 = 1 
for the total density, tthot = 0.3 for the hot component 
and fib = 0.013ft~ 2 for the baryonic fraction (Reeves 
1994). We also use ft = 0.5 for the dimensionless Hubble 
parameter and a = 0.78 for the rms fluctuation ampli- 
tude within a top- hat sphere of radius of 8ft -1 Mpc. In 
order to reduce the amount of shell crossing we filter the 
short-wavelength modes of the linear power spectrum 
using a Gaussian filter with radius R{ = 3.3ft" 1 Mpc. 
The model is normalized to COBE as in Gorski et al. 
(1994). We use this single model that has been shown to 
provide a good description of the cluster correlation and 
large-scale properties (cf. Klypin & Rhee 1994; Borgani 
et al. 1997). 

In TB95 it was shown that using large simulation 
volumes the variance in the dipole estimates is dom- 
inated by the observer-to-observer scatter within the 
same realization and not by the different realizations 
of the same dark matter model. Therefore, in what fol- 
lows, we will be using only a single realization of the 
model. 



3.3 Cluster Sample Extraction 

The simulation described above allow us to identify clus- 
ters with peaks of the dark matter density field and a 
subset of these peaks are used to generate a simulated 
cluster sample. In our previous analysis of cluster sam- 
ples (K96) , we assumed that there was a correlation be- 
tween the height of the density peak and the optically 
defined richness of a cluster. In this work, we assume a 
similar correlation between the height of the peak and 
the cluster X-ray luminosity. To simulate realistic flux- 
limited samples of X-ray clusters we need the number 
density of their parent population which can be con- 
verted into a threshold for peak selection. This is used 
to generate the "parent" volume-limited catalogue of 
clusters. Then by imposing a flux limit, and the shape 



5 

of the ROSAT selection function we extract realistic X- 
ray cluster samples. 

3.3.1 "Parent" Cluster Catalogue (C) 

By integrating the luminosity function from the proper 
lower to the upper luminosity limit, we can get an esti- 
mate of the mean number density of the "parent" dis- 
tribution of X-ray clusters: 

h c = J $(L) AL. (10) 

Taking the lower and the upper bound given in section 
3.1, we obtain n c « 3.3 x 10" 5 ft 3 Mpc" 3 . Thus the num- 
ber of "parent" clusters included within the cubic box is 
N c ss 28 900 and the number within an inscribed sphere 
of radius 480ft" 1 Mpc is ~ 15 100 which corresponds 
to a mean intercluster separation of d c « 31.3ft" 1 Mpc, 
which is characteristic of the APM type of clusters (Dal- 
ton et al. 1994). 

3.3.2 The mock 'ROSAT' sample (= R) 

To simulate a ROSAT-like flux-limited X-ray cluster 
sample we use the selection function, defined in equa- 
tion (^|), from which we can predict the correspond- 
ing number of X-ray clusters above the flux limit, 
Sum = 1.5 x 10~ 12 erg s" 1 cm -2 , lying within a shell 
between r and r + Ar, via: N(r) = iirr 2 n<f>(r)Ar. Fig. 1 
shows the resulting 4>(r) as a function of distance. No- 
tice the very long tail extending out to large distances. 
The maximum of N(r) turns out to be around z « 0.11, 
consistent with that predicted by the RASS (L. Guzzo, 
private communication). These considerations lead to a 
value of iV c w 1300 clusters within our limiting redshift 
of z ~ 0.16. We are therefore in perfect aggreement with 
the original survey target list of about 700-800 southern 
clusters (within z m 0.2). 

Finally we perform a consistency check, that of re- 
covering from our simulated flux-limited samples the 
observed N — S distribution (de Grandi 1996) which is 
well-described by a power law of the form 

N{> Si) ~ ASr a , (11) 

with q = 1.21 ± 0.15 and normalization constant 
A = 11.34+2 ^ sr" 1 (10~ n ergs" 1 cm" 2 ) . Using 500 
random observers, a flux-limit comparable to that of de 
Grandi (i.e. Sn m — 3 x 10" 12 ergs" 1 cm" 2 ) and fitting 
equation @, we find a = 1.22 ± 0.05, A = 11.7 ± 1 
in units as before and in excellent aggrement with de 
Grandi (1996). 

3.3.3 A 'test' sample (T) 

This test sample is constructed in order to study the ef- 
fects of missing nearby clusters which are too extended 
to be picked up by the SASS identification procedure 
(inability of the standard analysis software system to 
define them as such, simply because it is point source 
detection algorithm, therefore it can not deal properly 
with extended X-ray emission; but see also section 3.1 
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for additional flux correction). We use exactly the same 
requirements as for the R-sample, but now we adopt 
-Rmin = 7.5ft" 1 Mpc instead. The number of objects 
per simulation is only slightly affected, since we in- 
clude the few clusters that lie between 7.5ft _1 Mpc and 
45ft~ Mpc (according to their mean separation we ex- 
pect « 10 clusters within these limits). This particular 
lower limit is dictated by the procedure with which we 
have identified the observers; all local properties are de- 
fined by using a top-hat sphere of radius 7.5ft" Mpc 
(for more details see TB95) . Although such an inclusion 
increases insignificantly the total number of clusters, it 
does influence all dipole properties to a great extent, as 
expected from the 1/r 2 law of gravity and as we will 
quantify in our subsequent analysis. 



3.5 Weighting Schemes 

We finish this section by discussing the appropriate 
weighting schemes, to be used in equations (0) and 
(B), for the various samples. The C sample, being the 
"parent" population, is volume limited and thus does 
not suffer from selection effects, so we adopt the direct 
weight: Wi = r ; ~ 2 . For the R,T-samples, which are flux- 
limited, we use Wi w $ — 7";~ 2 to correct for the miss- 
ing objects. Note that this weighting scheme assumes 
that the unseen objects follow the same luminosity func- 
tion and that they are spatially, strongly correlated with 
the observed objects. Note that for the T-sample we 
use 4>(r) — 1 for distances between 7.5h~ Mpc and 
45ft" 1 Mpc. 



3.4 Observer selection 

From a starting list of 5000 random observers, identified 
on the grid, we will be selecting appropriate subsets with 
different characteristics depending on which particular 
issue we want to study. 

Firstly we select only those observers for which 
their cluster dipole amplitude, on large scales (~ 
180ft" 1 Mpc), exceeds their shot-noise dipole. This re- 
sults in 1500 such observers, which we call random ob- 
servers. We imposed this last criterion since we are inter- 
ested in studying dynamically 'active' regions, for which 
the dipole is well defined and a relative good diagnostic 
of the acceleration field. 

As a second subset we are interested in defining 
LG-like observers. Imposing the usual LG requirements 
(similar peculiar velocity and local density contrast; cf. 
Borgani et al. 1995; Moscardini et al. 1996) we generate 
an insufficient number of observers to allow a mean- 
ingful statistical analysis to be performed. A more use- 
ful constraint for studying the large-scale characteris- 
tics of the cluster distribution (which we also used in 
K96) is furnished by the the structure of the observed 
Abell/ACO cluster dipole, which is of a two-step form 
with a significant contribution from large scales (cf. 
Scaramella et al. 1991; Plionis & Valdarnini 1991; Bran- 
chini & Plionis 1996). In other words, we will just as- 
sume that our "parent" X-ray cluster distribution has 
similar dipole amplitude and shape as the Abell/ACO 
optical clusters. Branchini & Plionis (1996) starting 
form the redshift space cluster distribution they recov- 
ered cluster true distances and peculiar velocities via 
a dynamical algorithm (cf. Yahil et al 1991; Strauss et 
al 1992). Their mean reconstructed real space cluster 
dipole has the following features. 

(i) V c « 2400 ± 800 kms" 1 (2cr range) and 

(ii) the usual two-step shape, with a final convergence 
at ~ 180ft" 1 Mpc. 

This requirement reduces the number of available ob- 
servers in the simulation to 200, observers which we will 
hereafter call LG-like observers. 



4 ANALYSIS AND RESULTS 

We compute the moments of the cluster distribution for 
all 3 samples (C, R, T), in bins of 20ft" 1 Mpc width. We 
sum the contributions from all objects up to i? max = 
480ft" 1 Mpc using periodic boundary conditions. 



4.1 Monopoles 

We plot in fig. 2 the average random observer 
monopoles, for all three sets of clusters calculated ac- 
cording to equation (^), which as expected increase lin- 
early with r. For the flux-limited (R and T) samples, 
however, we observe a systematic early flattening at 
~ 100ft" 1 Mpc, which induces an overall underestima- 
tion of the "parent" population (C) monopole, although 
they are boosted up by using the appropriate correc- 
tion factor <j>~ 1 (ri). Their relative underestimation of 
the "parent" -population monopole, measured at depths 
ii > -Rconv, is S M « 10 and 20 per cent respectively for 
the T and R samples. This implies that the mean den- 
sity of X-ray clusters, calculated using the flux-limited 
sample monopole (via n — M/4irR), would be an under- 
estimate of their true underlying mean density, which 
would then provide an artificially lower estimate of /3 
in equation ((?]). We note that, up to 100ft" 1 Mpc, we 
expect on average ~ 50 — 60 clusters. 



4.2 Dipole as a function of distance 

To study the similarities of the different sample dipoles 
shape we use the usual correlation coefficient analysis 
(see equation 17 of K96). Furthermore to examine the 
accuracy with which the "parent" cluster distribution is 
traced by the flux-limited R and T samples, we mea- 
sure their corresponding velocity (dipole) fluctuations 
as follows. 



AVi,o(R) 



Vi(R)-Vc(R) 
Vc(R) 



(12) 



The subscript i corresponds to either R or T-samples 
and the subscript C corresponds always to our "par- 
ent" cluster sample (see also equation 18 of K96). Evi- 
dently, if the "parent" cluster population is adequately 
traced by flux-limited samples, we would expect to find 
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AVi,c(R) ~ V R, while in the case of missing contri- 
butions we would obtain AVi,c(^?) < 0. 

As we anticipated for the case of flux limited sam- 
ples (section 2.3), we have found that in many cases 
the raw dipoles (before correcting for the shot-noise ef- 
fects) diverge as a function of depth beyond the depths 
were their parent 3D cluster dipole has already con- 
verged. Therefore, we will present results based on both 
shot-noise correcting techniques discussed in section 2.3, 
although we expect that for some observers such a cor- 
rection will underestimate their true dipole amplitude. 



4-2.1 Random observers 

In most of the cases the "parent" sample (C) dipole 
amplitude is constantly higher than that of the T- 
sample. Note that, in both these cluster sets, the 'lo- 
cal region' up to 45ft~ 1 Mpc is included. The dipole 
shapes are quite similar up to -Rmax, which is evident 
from their relative shape correlation coefficient being 
~ 0.84 ± 0.13. However, the dipole amplitude of the 
C sample is underestimated by the flux-limited sam- 
ple, i.e. AVt,c(^) ~ —0.18 and -0.1 for the shot-noise 
corrected and uncorrected case resepectively. 

Comparing the C and R dipoles we find 
AV R ,c{R) ^ - 0.3 and -0.2 for the shot -noise corrected 
and uncorrected case resepectively, while their rela- 
tive dipole-shape correlation coefficient is 0.81 ± 0.12. 
This implies that in many occasions the C dipole is 
dominated by substantial 'local' contributions, within 
(50 — 60)ft -1 Mpc which the R dipole inevitably miss 
by construction and thus underestimates even further 
the "parent" population dipole. 

Examining the dipole amplitudes of the three clus- 
ter samples, we identify the following categories in de- 
creasing order of frequency. 

(i) Vc » Vt > Vb. in 60 per cent of the observers. 

(ii) Vc ~ Vt ~ Vr (within 10 per cent) in 30 per cent 
of the observers and 

(iii) Vr > Vc in 10 per cent of the cases. 

In fig. 3, we plot dipoles for the three samples as 
a function of radial distance. Panels (a) and (b) are 
representative of the first category, panel (d) of the sec- 
ond category and panel (c) is of the third case. Most of 
the R-dipoles underestimates the "parent" population 
dipole up to a distance of ~ 100ft _1 Mpc, where they 
actually start gaining amplitude. The final convergence 
comes at R CO nv ~ 200ft -1 Mpc. 

We also plot in fig. 4, the relative velocity fluctua- 
tions (AV). Note that, in the C versus T case the veloc- 
ity fluctuation, AVc,t{R), convergence to its final value 
at ~ 120ft -1 Mpc, whereas the situation is entirely dif- 
ferent for AVc,r(-R), as expected from the fact that the 
local region (< 45/i -1 Mpc) has been excluded in the R 
sample. Here at depths below 100ft -1 Mpc, we observe 
a very low value AVc,r(-R) ~ —0.6, which reaches its 
final value (w -0.3) at R « 180ft" 1 Mpc. 

It is important to note that, had we not corrected 
at all for the shot-noise effects (section 2.3) we would 
have found lower underestimates of the parent cluster 



dipole amplitude (i.e. AV(R) ~ —0.1 ~ —0.2) as ex- 
pected from the fact that our shot-noise correction is 
maximal (see discussion in section 2.3). In general, this 
indicates that the underestimation of the true underly- 
ing cluster dipole amplitude by the flux-limited samples 
is counteracted at some level by the artificial dipole am- 
plitude enhancement due to the shot-noise effects. The 
shot-noise correction is however very important in order 
to identify the true dipole convergence depth. 

4-2.2 LG-like observers 

Using equation (^) we find similar AV values as in 
the random observer case (i.e. AVr,c(-Rconv) ~ —0.16 
and AVr,c(-Rcoiw) ~ —0.3), which again represent an 
underestimate of the C-dipole. However, the correlation 
coefficients of the dipole shapes are systematically larger 
in this case than for the random observers. We find 0.9 ± 
0.1 for the C versus T samples and 0.87 ± 0.09 for the 
C versus R ones. It is clear that, the three dipoles are in 
good agreement, as far as their shapes are concerned. 

We can break down the behavior of the dipole am- 
plitude for the LG-like observers into the same cate- 
gories as for the random observers case. Only slight 
changes occur in the relative frequencies with most 
notable the slight increase of the observers for which 
their flux-limited samples provide good estimates of the 
"parent" -population cluster dipole (C) (~ 40 per cent 
of total). We do not present plots for the LG-like ob- 
servers due to their qualitative similarity with those of 
the random observers (figs 3 and 4). 

Finally, we report that for both sets of observers we 
have also applied the classical shot-noise correction of 
section 2.3. Although we do not show any plots for rea- 
sons of brevity, results clearly suggest that whatever the 
adopted method for such corrections, both flux-limited 
cluster samples (R, T) underestimate the "parent" (3D) 
dipole by at least 10 per cent independently of the ob- 
server. 



4.3 Dipole alignment 

In fig. 5 we present the relative dipole misalignment 
angles, A9, between the three sample dipoles, averaged 
over the ensemble of all observers. We define A#c,t to 
be the misalignment angle between the dipole vectors 
of the C and T samples, and A#c,r between C and R. 
Panel (a) shows the A9 as a function of distance for 
the random observers while panel (b) for the LG-like 
observers. 

For the random observers, we obtain mean values of 
A# C ,T ~ 18° ± 15° and A0 c ,r ~ 47° ± 30°. The above 
numbers have been measured at fl conT . The big differ- 
ence of « 30° in the mean misalignment angles between 
A#c,t an d A#c,r can be attributed to the fact that the 
R dipole ignores (by construction) the local region (< 
45ft - Mpc). This again shows how important the 'lo- 
cal region' actually is in determining the cluster dipole. 
From panel (a) of fig. 5, it is clear that the misalignment 
angle between the C and T sample dipoles converge to 
its final value at ~ 100ft -1 Mpc while between the C 
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and R-dipoles this scale is ~ 200/i 1 Mpc. For distances 
well below this limit, we get AO ~ (60° - 70°). 

As seen from panel (b), the situation is greatly im- 
proved, when we use the LG-like observers, since by 
construction their dipole properties are not dominated 
only by the 'local region'. The large misalignment angles 
are reduced rapidly and converge to a significantly lower 
value at ~ 150/i _1 Mpc. We obtain A6> c ,t » 17° ± 11° 
with a median of 18° and A6» c ,r « 30° ± 20° with a 
corresponding median of 31°. 

Note, finally, that we did not apply any correction 
for shot-noise effects in the dipole direction analysis. A 
shot-noise correction is expected to reduce the dipole 
misalignment angle by at least 10°. 

4.4 A Cluster Bias Estimate 

Using equation (Q) and the input value of fi (= 1), 
given by the model we simulate, we can estimate the 
cluster bias parameter, via: 



6c 



1\ 

-y 



\v\cAR) 
\u\i(R) ' 



(13) 



where N is the total number of observers used and 
R > -Rconv ■ The small subscript c indicates cluster veloc- 
ity for the three samples (see section 3.3 and equation 
7), whereas i shows the specific observer used in the 
summation. A more elaborate version of this procedure 
was tested and found to provide a reliable estimate of 6 C 
in K96. Using equation (w|) and those LG-observers for 
which their peculiar velocity vector lie within ~ 25° of 
their measured cluster dipole we obtain for the "parent" 
sample: 



b c = 4.8 ± 1.5 , 



(14) 



a value which is very similar to that derived by the ratio 
of the cluster to dark-matter variances (measured on 
scales > 50/t -1 Mpc). However the method of equation 
( [l3| ) produces a large observer-to-observer scatter (6 to 
8 times larger than in variance method) which is the 
outcome of the fact that the very local contributions to 
the dynamics of some observers may not be represented 
in the cluster distribution, since these can be due to the 
very local inhomogeneities in the matter distribution. 

The underestimation of the "parent" -population 
dipole amplitude by the flux-limited samples will in- 
evitably reflect itself on an artificially lower bias pa- 
rameter, given by equation (|l3|). In fact we obtain 
6 C « 3.95 ± 1.6 and 3.45 ± 1.4 for the T and R sam- 
ples respectively. 

4.5 Summary of Moment Analysis 

The flux-limited cluster samples underestimate the X- 
ray "parent" population dipole. All dipoles are in good 
agreement (within 10 per cent), in about only one in 
three cases. The main effect of imposing a flux-limit 
is a general loss of dipole amplitude, which amounts to 
10/20 per cent for the T sample and up to 20/30 per cent 
for the R sample (depending on whether one corrects or 
not for shot-noise effects). Note that the above results 



are independent of the two methods used (discussed in 
section 2.3) to correct for the shot-noise effects. The 
relative difference between the two flux-limited sam- 
ples (which can reach a maximum amplitude of ~ 15 
per cent) is due to the exclusion of the 'local region' 
from the R sample, which corresponds to the inability 
of the RASS to sample the local extended X-ray clus- 
ter sources. The misalignments between the flux-limited 
sample dipoles and that of the "parent" -population are 
quite large, ranging from 20° to 47° on average (for the 
T and R samples respectively). 

If we restrict our analysis to those observers that 
have a "parent" cluster population that exibits a similar 
dipole to the observed Abell/ACO cluster case (i.e. hav- 
ing significant contribution from large scales; cf. Bran- 
chini & Plionis 1996), the two flux-limited samples have 
dipole directions in much better aggreement with that 
of the "parent" cluster dipole (AO ~ 17° and 30° respec- 
tively for the T and R samples). The dipole amplitude 
underestimation by the flux-limited samples is however 
of a similar amplitude as in the random observers case. 

Note, importantly, that the underestimation of the 
"parent" cluster dipole amplitude by these flux-limited 
samples compared would lead, using equation (Q) , to an 
incorrectly higher estimate of the f3 c parameter and thus 
to either an artificially higher value of fi for a specific 
value of 6 C or to an artificially lower value of 6 C for a 
specific value of fi Q . 



5 X-RAY CLUSTER CONTRIBUTION TO 
XRB 

5.1 A brief overview 

Since it was brought forward by Giacconi et al. (1962), 
the XRB has been something of an enigma for theo- 
retical cosmology. Its origin and its major constituents 
are still debatable nowadays, and it is still unclear how 
much of it is contributed by resolved sources and how 
much by a discrete component (see review by Zamorani 
1993 as well as Hasinger et al. 1993; Hasinger 1992; 
Hasinger 1995). The diffuse hypothesis (thermal breem- 
stralung hot plasma with T ~ 40 keV) has largely been 
abandoned, because such a hot component would gen- 
erate severe distortion on the CMB spectrum, which 
is not detected by COBE. The alternative solution of 
assuming that the XRB is made up of the integrated 
X-ray flux from discrete objects is more attractive, but 
as yet no specific class of sources has been found with a 
spectrum identical to that of the XRB. 

Setting the spectral difficulties to one side, we shall 
attempt to quantify the contribution to the XRB from 
clusters of galaxies. Many similar attempts have been 
made for various other sources (active galaxies, QSOs) 
as well as clusters, in both soft (E < 2 keV) and hard 
(E > 2 keV) energy bands; see Shanks et al. (1991) 
(and references therein), Blanchard et al. (1992), Geor- 
gantopoulos et al. (1993), Boyle et al. (1994), Roche et 
al. (1995) (and references therein) for the QSO contri- 
bution to the soft XRB, Lahav et al. (1993) (and ref- 
erences therein) discuss the contribution from galaxies, 
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Anvi (1978) and Maccacaro et al. (1984) have discussed 
the case for AGN and Piccinotti et al. (1982), Gioia et al. 
(1990), McKee et al. (1980), Schwartz (1978), Henry et 
al. (1992), Soltan et al. (1995), Kitayama & Suto (1996), 
Kitayama, Sasaki & Suto (1997) and Oukbir, Bartlett 
& Blanchard (1997) are mostly concerned about X-ray 
clusters. 

We will attempt to estimate the contribution of X- 
ray clusters, having a X-ray luminosity function given 
by equation (^) to the soft XRB using two simple 
methods (see Maccacarro et al. 1984, their section 4), 
not worrying for the different cosmological backgrounds 
that could affect the detailed modelling of such con- 
tribution (cf. Kitayama & Suto 1996; Kitayama et al. 
1997). Note that the present application holds only for 
a flat model with vanishing A term, since this is the 
background that de Grandi (1996) has assumed in de- 
riving the relative parameters of the ESO-KP X-ray 
luminosity function. 



5.2 Luminosity Function Method 

First we compute the contribution of X-ray clusters ei- 
ther in terms of flux per steradian or in terms of lu- 
minosity density (ergs -1 frMpc -3 ). In other words, we 
compare either intensities of clusters and XRB, or vol- 
ume emissivities at a given energy band. The contri- 
bution of clusters is estimated integrating the product 
L$(L): 



Jc 



L$(L) dL, 



(15) 



where L m i n ~ 4 x 10 41 h~ 2 ergs -1 and L max ~ 
10 45 h~ 2 erg s -1 (see section 2.2). The quantity J c ex- 
presses the integrated cluster volume emissivity. Since 
there is no clear picture which describes the correspond- 
ing volume emissivity of the XRB in such a soft band, 
we will simply be using the well-defined intensity of 
the XRB in the 1-2 keV band, namely 1% ~ 1.25 x 
10 -8 ergs -1 cm -2 sr -1 of Hasinger et al. (1993), as the 
basis of our estimate. Assuming values for the spectral 
index of the soft XRB in the usual range of 1-1.2, we first 
convert the latter value in the monochromatic band (at 
1.5 keV) and then we use a conversion factor to trans- 
form it to the 0.1-2.4 keV one. Such a calculation yields 
Ix « 4 ± 0.2 x 10 -8 ergs -1 cm -2 sr -1 (0.1-2.4 keV), 
for the XRB intensity in the usual ROSAT band. Evi- 
dently, we should compare intensities with luminosities 
in identical bands and for this purpose we either convert 
the cluster volume emissivity in flux, or the intensity of 
the XRB in light density units. The former requires an 
accurate knowledge of the spectral index of X-ray clus- 
ters, for the conversion to be valid. We will assume a 
typical value for the spectral index of X-ray clusters in 
the soft band of the order of 0.4 (photon index 1.4). 
The relation between volume emissivities and intensi- 
ties (Schwartz & Gorsky 1974; Anvi 1978; Soltan et al. 
1995) is given by 



The subscript i corresponds to the objects under study; 
one can use equation (16) interchangeably for h or Ji. 
The look-back integral depends on parameters such as 
H , q and Oi (a zero cosmological constant is assumed) . 
Varying q a in the range between 0.1 and 0.5 does not af- 
fect the the look-back time integral by more than 10 per 
cent. The cluster contribution to the soft XRB, if q a is in 
the above range, is therefore well within the uncertain- 
ties quoted throughout this section. Furthermore the de- 
pendence on Ho drops out from the final result. Using 
the spectral index of the XRB given above, we can solve 
for Jx. The upper limit of the integral presents the max- 
imum distance up to which, sources generate the XRB. 
Values between 3, 4 and oo do not influence the final re- 
sults by more than 3 per cent. In what follows we will be 
using z max = 4 and q = 0.5. We then obtain from equa- 
tion (16) that J x « 4.1 ± 0.1 x 10 39 ergs -1 /iMpc -3 
(0.1-2.4 keV). Adopting the corresponding value for the 
clusters (equation 15), we can express the fractional con- 
tribution of X-ray clusters to the total of the soft XRB 
as 

Co.l-2.4kcV = ^ ~ (9 ±1)%, (17) 
Jx 

which is independent of the value of H - 

We can obtain another estimate of Co.i-2.4keV, 
starting with the formalism of Schwartz & Gorsky 
(1974), using equation (16) and the spectral intensity 
of the XRB for energies between 1 keV and 21 keV and 
then solving for Jx. The spectral intensity is expressed 



1(E) « 8.5 E- 



keV 



(18) 



Ii = 



c Jj 
AttHo 



(1+*) 



(l + z) 3 Vl + 2qoZ 



:dz. 



(16) 



ke V cm 2 s sr ' 

for 1 < E < 21 in keV's. It is apparent though, that 
such a scheme does not cover the lower limit of our soft 
band (at 0.1 keV). We can circumvent, however, this 
difficulty by using the following alternatives. The first 
is to extrapolate the shape of the spectrum at low en- 
ergies, so as to match our lower bound, and the second 
is to perform calculations in the perfectly allowed (1-2 
keV) band, which we later convert in the soft ROSAT 
one as before. We choose the second method, since it 
does not violate the energy range and it appears to be 
more realistic. However, the former method produces 
figures always within ± la of the latter. We caution the 
reader that there is no well-known 1(E) for the band 
we are interested in and this energy range is mostly un- 
explored in terms of spectrum. Taking again the same 
set of parameters as before (H , q , ax, Jx), we obtain 
Jx ~ 3.95 ± 0.1 x 10 39 ergs -1 /iMpc -3 (0.1-2.4 keV), 
which is in good agreement with the value previously 
computed and which gives Co.i-2.4kcV ~ (10 ± 1) % . 
Note that the errors of both estimates reflect the uncer- 
tainty of the choice of the spectral index of the soft 
XRB, being used in equation (16) (cf. Soltan et al. 
1995). 

A further estimate of Co. 1-2.4 k C v is obtained us- 
ing an estimate of the spectral index of clusters of 
the order of a c — (0.4 — 0.5), as in Maccacaro et 
al. (1988), and using the value of J c from equation 
([[5]) in the 0.1-2.4 keV band, we can then solve equa- 
tion (|l6|), this time for the total intensity for X-ray 
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clusters, I c , in the same energy range. We find that 
J c « 4.5 ± 0.1 x lO^ergs^cm^sr -1 (0.1-2.4 keV). 
By comparing this with the previous estimate of 7x (0.1- 
2.4 keV), we find 

Co.i-2.4k.v = -r- « (11 ± l) 96 ( 19 ) 

iX 

consistent with the previous estimates. For all the above 
estimates, we have assumed that the spectral indices of 
both clusters and XRB are not subject to any change 
when passing from one energy band to another and al- 
though a few such assumptions have been used, this 
framework does yield consistent results. 

5.3 Number-Flux Distribution Method 

A different approach is used here, which is independent 
of the value of q a , that enters in the evaluation of the 
look-back integral. We assign a random luminosity to 
each simulated cluster using the X-ray luminosity func- 
tion of equation (^) , which we then convert to a flux us- 
ing its distance. For this, we use the constant multiplier 
of de Grandi (1996) to transform 0.5-2 keV fluxes (or lu- 
minosities) to the broad 0.1-2.4 keV band (as discussed 
in section 2.2). Note that the particular CHDM simula- 
tion we use is irrelevant to this analysis; we could have 
even used a Monte-Carlo cluster distribution provided 
that the N(r) distribution, predicted by the X-ray clus- 
ter selection function, would be preserved. We plot in 
fig. 6 the distribution of X-ray luminosities as a function 
of distance. The solid line corresponds to the imposed 
flux limit of 1.5 x 10 -12 ergs -1 cm -2 , for any given dis- 
tance. The plot shows results for a single observer; we 
have verified that this behavior of the distribution is 
entirely independent of the observer position. 

Having obtained the fluxes of the ~ 1300 X-ray- 
like clusters, we simply add their individual contribu- 
tions towards the total cluster X-ray emission in the 
corrected soft band. To estimate the scatter due to 
sampling effects we use 500 random observers from the 
initial ensemble and we repeat the procedure for each 
number-flux relation. We report results, as usual, as 
the mean and the ± la errors. We find that, 7 C ~ 
4.25 ± 0.3 x 10~ 9 ergs- 1 cm- 2 sr- 1 (0.1-2.4 keV). We 
compare intensities of clusters and XRB as computed in 
section 5.2. We obtain Co.i-2.4kev ~ 10.5 ± 1 per cent 
for the contribution of the discrete set of these clusters. 
Alternatively, we can again assume equation ([Hi]) and 
the same spectral index for clusters, in order to con- 
vert fluxes per unit area to volume emissivities. This 
approach yield identical results. 

However, since we have used our simulated cata- 
logue which is limited to a radius of 480/i -1 Mpc, we 
need to take into account the contribution of the extra 
~ 300-400 clusters predicted by the X-ray cluster selec- 
tion function to lie z ^ 0.16, which we estimate to be no 
more than 1.5 per cent of the previously quoted value. 
Therefore this estimate is in very good agreement with 
our previous findings (see equations [l?] and |l^) . 

We, therefore, quote the contribution of clusters to 
the soft XRB, as an average of the two methods used 
above. This yields 



Co.i-2.4kov « (10 ±2)%, (20) 

with quoted uncertainty covering the la deviations of 
all previously computed averages. A similar value was 
obtained recently by Oukbir et al. (1997). 



6 CONCLUSIONS 

We have used numerical simulations of the cluster distri- 
bution and the de Grandi (1996) X-ray cluster luminos- 
ity function, to create samples of X-ray type clusters. 
We investigate their relative dipole properties, by pay- 
ing particular attention to what extent the underlying 
"parent" population cluster dipole can be recovered by 
the flux-limited subsamples. 

We analysed the results based on 2 different sets of 
observers; random observers and those that have a clus- 
ter dipole structure in reasonable agreement with that of 
the Abell/ACO clusters (LG-like observers). We found 
that on average the "parent" cluster dipole is underes- 
timated by ~ 20/30 per cent in the former type of ob- 
servers and slightly less, ~ 10/20 per cent, in the latter 
type (for the raw and shot-noise corrected cases respec- 
tively). However, the profile of the "parent" dipole am- 
plitude as well as the dipole direction is much better re- 
covered in the LG-like observer case (AO ~ 30° and 47° 
respectively for the LG-like and random observers). We 
also found that the exclusion from the flux-limited X- 
ray cluster samples of the local region (^ 45/i _1 Mpc), 
as is the case with RASS, will result into larger under- 
estimates of the "parent" dipole amplitude and larger 
misalignments of the dipole directions. 

The generic loss of dipole amplitude will be im- 
portant in interpreting cosmological information coming 
from samples of X-ray selected clusters, such as their 
relative bias factor and the inferred values of Qq. 

Finally, we attempted to estimate the contribution 
of these X-ray clusters to the soft (0.1-2.4 keV) XRB for 
the particular case of a flat universe with Qa ~ 0. We 
use the standard method for deriving the cluster volume 
emissivity, as well as that using the cluster number-flux 
relation. The fractional contribution of clusters to the 
soft ROSAT energy band is estimated to be 10 ± 2 per 
cent, using both techniques. 



Note added in Manuscript: When this project 
was in its final stages, two new X-ray flux-limited clus- 
ter samples appeared in the literature, namely the X- 
ray Brightest Abell-type Cluster sample (XBACs; Ebel- 
ing et al. 1996b) and Brightest Cluster Sample (BCS; 
Ebeling et al. 1997a and Ebeling et al. 1997b). The 
XBACs sample consist of 242 Abell/ACO X-ray clus- 
ters above the limiting flux of 5 x 10 -12 erg s _1 cm" 2 
and it is the biggest all-sky X-ray cluster list, up till 
now. The BCS (largest in the northern hemisphere) 
comprises 199 X-ray Abell, Zwicky as well as other 
purely X-ray selected systems, being flux-limited above 
4.45 x 10 -12 ergs -1 cm -2 . Both samples are selected at 
high galactic latitudes ([ b\ £ 20°) and in the same soft 
energy band (0.1-2.4 keV). In both samples the cluster 
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identification procedure is mostly based on the Voronoi 
Tesselation Percolation (VTP) algorithm (see also Allen 
et al. 1992; Ebeling & Wiedenman 1993; Crawford et al. 
1995), which not only samples nearby extended X-ray 
emission much more adequately than the SASS algo- 
rithm, but also increases significantly the accuracy with 
which fluxes are computed. This, together with the fact 
that their depth coverage is large (z ~ 0.2 and 0.3 for 
XBACs and BCS respectively), renders both samples 
ideal tools for large scale structure studies. These cata- 
logues together with the forthcoming ESO KP will def- 
initely constitute three of the more significant X-ray 
databases for future cosmological research. We, there- 
fore, intend to study these samples extensively along 
the lines set in this work. 
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FIGURE CAPTIONS 



Figure 1. X-ray selection function for a model 
ROSAT-like cluster population. 

Figure 2. The average, over 1500 random ob- 
servers, monopole as a function of radial distance for the 
three samples of clusters. Errorbars reflect the observer- 
to-observer scatter. Filled dots are for "parent" (C) 
monopoles, open dots are for T monopoles and filled 
triangles refer to R monopoles. 

Figure 3. Dipole amplitude as a function of radial 
depth for some characteristic random observers. The 
solid line corresponds to the R, open circles to the T 
and filled dots to the "parent" (C) samples. The up- 
per 2 panels represent the most frequent cases, lower 
right those less frequent and finally the lower left are 
the rarest of all. 

Figure 4. Relative velocity fluctuations (equation 
14) for the 1500 random observers. Errorbars reflect the 
observer-to-observer la scatter and are shown only for 
one set of points for clarity. Filled dots correspond to 
fluctuations between the "parent" and T samples while 
open dots between the "parent" and R samples. 

Figure 5. Relative misalignment angles, A8. ran- 
dom observer results are shown in panel (a) while those 
of the LG-like observers are shown in panel (b). Error- 
bars are only shown for one set of points, in each panel, 
for clarity. In both panels filled symbols correspond to 
angular fluctuations between C and R samples, whilst 
open symbols represent misalignment angles between C 
and T populations. 

Figure 6. Luminosity-distance diagram for a 
model ROSAT-like sample, as viewed by a typical ob- 
server. The solid line corresponds to the cut-off intro- 
duced due to the flux limit imposed. The plot shows the 
X-ray cluster distribution of ~ 1300 such objects seen 
by a specific observer, but the result (shape) is indepen- 
dent of the observer's choice. 
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X-ray selection function 




Wonopole amplitude 



_ Filled circles: C — sample 
Open circles: T— sample 
Filled triangles: R — sample 
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